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Abstract 

The influence of intense magnetic fields on the behavior of color superconducting cold quark 
matter is investigated using an SU(2)j NJL-type model for which a novel regulation scheme is 
introdnced. In such a scheme the contributions which are explicitly dependent on the magnetic 
field tnrn out to be finite and, thns, do not require to be regularized. As a result of this, non¬ 
physical oscillations that might arise in the alternative regnlarization schemes previonsly nsed in 
the literature are naturally removed. In this way, a clearer interpretation of the physical oscillations 
is possible. The sensitivity of our results to the model parametrization is analyzed. 
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I. INTRODUCTION 


At asymptotically large chemical potentials, the fact that cold quark matter behaves as 
a color superconductor can be shown by using perturbative methods in the context of quan¬ 
tum chromodynamics (QCD)[l|. However, such methods cannot be applied in the range 
of moderate densities relevant for, amongst others, the astrophysics of strongly magnetized 
compact stellar objects known as magnetars. We recall here that, although it is generally 
accepted that these objects can have surface magnetic helds up to 10^^ G|^, the estimates 
for the magnetic held values at their centers are model dependent to some extent, ranging 
between B ~ 10^® —10^° G (see e.g. Refs. 3-^). Since the well-known sign problem prevents 
lattice QGD calculations from being performed at sufficiently low temperatures and hnite 
chemical potential, one has to rely on effective models to analyze the behavior of magnetized 
quark matter in this region. One particular model that has been extensively used for this 
purpose is the Nambu-Jona-Lasinio (NJL) model This is an effective model originally 
devised to study the dynamics of chiral symmetry breaking, in which gluon degrees of free¬ 
dom are integrated out and interactions are described by local four-quark interactions. The 
incorporation of additional di^ark interactions into the model allows for the description of 
color superconducting matter [7|. In this context, the effect of a constant magnetic held has 
been analyzed by several authors 8l-ll3l|. At this point, it is important to remark that the 
local character of the interactions considered in the NJL-type models leads to divergences in 
the momentum integrals which need to be handled in some way in order to completely dehne 
the model and yield meaningful quantities. Several regularization procedures are possible 
even in the absence of magnetic helds Moreover, when magnetic held is introduced, the 
vacuum energy acquires a Landau level (LL) structure and an additional care is required in 


the treatment o: 
reported in Ref. 


t 


le divergences. An elegant way of treating the regularization has been 


14j | for the model in the absence of color superconductivity. The procedure 


follows the steps of the dimensional regularization prescription of QGD, performing a sum 
over all Landau levels in the vacuum term. This allows to isolate the divergence into a term 
that has the form of the zero magnetic held vacuum energy and that can be regularized in 
the standard fashion. It should be stressed that similar expressions for the magnetic held 
dependent terms can be obtained using a method based on the proper-time formulation 151. 


So far, however, this procedure has not been applied to the case in which color superconduc- 
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tivity is present. Instead, existing calculations j9l-ll3| remove the divergences by introducing 
some type of regulator function for each Landau level separately. This procedure, however, 
might in general introduce unphysical oscillations. A discussion on this can be found in 


Refs. 


16 


18| . where it is also observed that the use of smooth regulator functions improve 


the situation. In fact, this allows to identify possible physical oscillations appearing in some 
cases { 9 , 101. However, an even clearer interpretation of the results could be obtained if the 
unphysical oscillations were removed altogether with another scheme, especially at hnite 
chemical potential and in the presence of color superconductivity. The main purpose of this 
work is to investigate the influence of a constant magnetic held on cold superconducting 
quark matter in the framework of the NJL-type model, using a regularization procedure in 
which the contributions that are explicitly dependent on the magnetic held turn out to be 
hnite and, thus, do not required to be regularized. This procedure will be referred to as 
“Magnetic Field Independent Regularization” (MFIR), and it can be considered an exten¬ 
sion of the method described in e.g. Ref.jl^ to the case in which color pairing interactions 
are present. Since the aforementioned unphysical oscillations are completely removed in this 
scheme, we will complement our analysis by performing a detailed study of the resulting cold 
matter phase diagrams, including their dependence on the parameters of the model and, in 
particular, the coupling strength of the diquark interactions. Being mostly concerned with 
the issues related to the model regularization procedure we will, for simplicity, assume all 
the quark species to have a common chemical potential leaving the incorporation of the 
neutrality and (3 equilibrium conditions relevant for stellar matter applications for future 
studies. 


We organize the article as follows. In Section II we present the Nambu Jona Lasinio 
model with magnetic held and diquark interactions. In particular, we briehy review the 
regularization schemes used in the literature and describe in some detail the MFIR scheme 
introduced in this work. The model parameters used in our numerical calculations are also 
given. In Sect. Ill we compare our results for the behavior of the cold and dense magnetized 
quark matter with those previously reported in the literature. In Sect. IV we present results 
for the phase diagrams in the eB — fi plane as obtained using MFIR for different interaction 
coupling ratios and parameter sets. In Sect. V we present our conclusions. Finally, in the 
Appendix, several details of the formalism of the MFIR scheme are described. 


3 









II. MAGNETIZED COLD QUARK MATTER WITHIN THE SU(2)/ NJL MODEL 
IN THE PRESENCE OF COLOR PAIRING INTERACTIONS 


A. The thermodynamical potential in the mean field approximation 


+6- 


We consider a NJL-type SU(2)/ Lagrangian density which inclndes scalar-psendoscalar 
and color pairing interactions. In the presence of an external magnetic field and chemical 
potential it reads: 

C = Ip i p — rric + fJ. pJ 

[Ppp)f + + H [{i^p^ efel'^^ip){iip ejel'^^ip^)] . (1) 

Here, G and H are conpling constants, tp = (n, d)^ represents a qnark field with two flavors, 
xp^ = Ctp'^ and tp^ = ip'^C, with C = * 7 ^ 7 °, are charge-conjngate spinors and r = (ri, T 2 , T 3 ) 
are Panli matrices. Moreover, = (ec)^“^ and are antisymmetric matrices in color 

and flavor space respectively. Fnrthermore, rric is the (cnrrent) qnark mass that we take 
to be the same for both flavors and /i is the qnark chemical potential. The conpling of 
the quarks to the electromagnetic field Afj, is implemented through the covariant derivative 
— ieQA^. Note that here we are dealing with “rotated” fields. In fact, as is well 
known, in the presence of a non-vanishing superconducting gap A, the photon acquires a 
finite mass. However, as shown in Ref. |l9 1. there is a linear combination of the photon and 
the eighth component of the gluon field that leads to a massless rotated U{1) field. The 
associated rotated charge matrix Q is given by 

As 


Q — Qf 


2x/3 


( 2 ) 


where Qf = diag(2/3, —1/3) and is the color quark matrix Ag = diag(l, 1, —2)/\/3. 
Then, in a six dimensional flavor-color representation {ur, Ug, u^, dr, dg, dh), the rotated q for 
different quarks are: Ur = 1/2, Ug = 1/2, Ub = l,dr = —1/2, dg = —1/2, db = 0. The rotated 
unit charge e is given by e = ecos 6 , where 6 is the mixing angle which is estimated to be 
~ 1/20 20|. In the present work we consider a static and constant magnetic field in the 
3-direction, Ag, = 6 g 2 XiB, which in fact is also a mixture of the electromagnetic field and 
color fields. 

In what follows we work in the mean field approximation (MFA), assuming that the only 
non-vanishing expectation values are < ipip >= —(M — mc)/ 2 G and < i'lp^ efe/.'jb'P’ >= 
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—A/2H, which can be chosen to be real. Here, M and A are the so-called dressed qnark 
mass and snpercondncting gap, respectively. The resnlting MFA thermodynamic potential 
at vanishing temperatnre reads 


f^MFA — 


(M 


rur 


AG 


-w- ^ 

kl=o,^,i 


(3) 


where 

P\q\=o 

P\q\ = l 

P\q\ = l/2 

Here, we have introdnced ak = 



2 — 6 ko and 


Et 


\/pI + 2keB -|- ± jj, 


\/pI + keB + ± /i -1- A^. 


(4) 

(5) 

( 6 ) 


(7) 


Clearly, Eqs.dMH]) are divergent and, thns, reqnire to be regnlarized. Some alternative 
schemes to achieve this will be discnssed in the following snbsection. Given the corresponding 
regnlarized form the associated gap eqnations for M and A then read 


‘MFA _ Q 

d{M, A) 


( 8 ) 


For each valne of /i and eB, several solutions of these equations will generally exist, corre¬ 
sponding to different possible phases, and the most stable solution is that associated to the 
absolute minimum of the thermodynamic potential. 


B. Regularization schemes 


As already mentioned, the contributions from Eqs.(jl}lH]) need to be regularized. In pre¬ 


vious studies jo-ll3| this was accomplished by introducing some cutoff function hA{q) in the 
corresponding integrands, with q = p in the case of Eq.(jl]) and q = = \/pI + 2k\q\B 

for Eqs. fl5]l6|) . The obvious and simplest choice would be to take h^{q) = Q{A — q). In 
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what follows, we will refer to this regularization scheme as “sharp function regularization” 
(ShFR). In this case, the integral in Eq.(jl]), which is magnetic held independent, is cut off 
when p = A. The contributions coming from Eqs.(|MS]) include a sum over Landau levels, 
and the integral in each of these is cut off for the that satishes A = , that is, 

such that the momentum and magnetic held contribution to the quark dispersion relation 
does not exceed the value of the cut-oh. This would seem like a natural way of extending 
the 3D sharp cutoh zero magnetic held regularization to the hnite eB case. However, the 
magnetic held dependence of the prescription brings in strong unphysical oscillations. To 
minimize the ehects of this magnetic held dependence, the aforementioned studies have been 
carried out replacing the Heaviside function with a smooth regulator. For dehniteness in 
this work, we consider the function h\{q) = 1/(1 + exp[(g/A — l)/a]). We have verihed that 
other possible choices lead to similar results. To choose the value of the constant a that 
determines the regulator smoothness one is limited by the fact that a too steep function does 
not improve over the ShFR results and that a too smooth function leads to values of the 
quark condensate in absence of the mag netic held which are quite above the phenomenolog¬ 
ical range. Here, we follow Refs. jlT. 12 | and consider a = 0.05. In what follows we will refer 
to this regularization scheme as “smooth function regularization” (SmFR). The similarities 
and diherences in the results yielded by the two regularization schemes introduced so far 
will be discussed in the following section. There it will be noted that, although somewhat 
suppressed, the undesired oscillations are still present in SmFR. It is important to remark 
at this stage that P|q|=o,i can be rewritten in terms of a vacuum and a matter contribution, 
of which only the hrst one is divergent. Thus, there is certain ambiguity on whether the 
regularization function has to be included in the matter term or not. Having analyzed both 
possibilities, we verihed that in our case, this ambiguity amounts at most to small quantita¬ 
tive differences. The results to be presented correspond to the case where only the vacuum 
energy is regularized. 

To fully get rid of the above mentioned regularization artifacts, we introduce in what 
follows an alternative scheme in which the contributions that are explicitly dependent on 
the magnetic held turn out to be hnite and thus do not need to be regularized. We will refer 
to this regularization scheme as “magnetic held independent regularization” (MFIR). We 
start by considering the P\q\=o contribution. Since it is independent of the magnetic held, it 


can be treated in the usual way 


Introducing a sharp 3D cutoh we get 
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( 9 ) 


Piii\=o = —, I <ip p‘ 

VT ./n 


0(/i - M) 


TT^ 


/i(/i2 - M2)3/2 (^2 _ 


-/l 


M 


v//i2 - M2 


where h(z) = (2 + + ;24in[^/Vl + z2]. 


In the case of P\q\=i we note that, except for the specihc value of the quark charge, the 
con.espo„di„g explion co.c.des with that aaaly.ed Ref.fl where no color pairing 


interactions were considered. Following the steps discussed in that reference we get 

^max _ _ / I t O 

+ V 

y (VI. — .s, — .St III I — 

47r2 


TT 


+ 




P\q\ = l = — dpp 


Oik 


>0 


k=0 


p\ p 


4 In 


Sk 


{ePy 


27r2 


2 2 ' 

n-l,a) + ^^ lna + ^ 


( 10 ) 


where x = M‘^/{2eB), kmax = Floor[(/i^ — M'^)/{2eB)] and Sk = + 2keB. In Eqs.dU]- 

[IQD, the first term is a vacuum contribution which does not explicitly depend on the magnetic 
held and the second term is the matter contribution. The last term in Eq. ffTOj) is the explicit 
magnetic held contribution to the vacuum, which has been isolated into a hnite term. 

The case of |g| = 1/2 is more involved. However, as discussed in detail in the Appendix, 
it can be cast into the form 


-^ 51 = 1/2 dpp^ + E^) + 

7i Jn 


27r2 


cu 1 \, y~y^^ I 2 /^ 


{iBy 


27r2 


dp 


OO ^(DO 

/(4 + fc) - 2 / dxf{p‘^ + x) 
k=0 do 


( 11 ) 


where {\Jp'^ + M2 ± /i)2 + A2, y = (M^ + A^)/(eH) and 

/(4 = X] \J {•\/z + 2x + s p/ViBy + y-2x--s/z + y . (12) 

s=±i L 

In this expression a 3D sharp cutoh has been introduced to regularize the hrst term, i.e. the 
one that contains contributions from the vacuum and matter which do not explicitly depend 
on the magnetic held. Note that, as in the case of vanishing eB discussed in e.g. Ref. 21| . 
these cannot be disentangled into two terms unless A = 0. The second term is the vacuum 
magnetic contribution analogous to the |g| = 1 case. Finally, the third term is an additional 
explicitly magnetic held dependent matter contribution which, as shown in the Appendix, 
turns out to be hnite. 
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C. Model parametrization 


In order to analyze the dependence of the results on the model parameters, we will 
consider two SU(2)/ NJL model parameterizations. Set 1 corresponds to that leading to 
Mq = 340 MeV while set 2 to that leading to Mq = 400 MeV, within the MFIR regular¬ 
ization. Here, Mq represents the vacuum quark effective mass in the absence of external 
magnetic helds. The corresponding model parameters are listed in Table [H 

TABLE I. Parameter sets for the SU(2)j NJL model. In both cases, empirical values in vacuum 
for the pion observables are reproduced, ttItt = 138MeV and /tt = 92.4MeV. 


Parameter set 

Mq 

me 

GM 

A 

— < uu 


MeV 

MeV 


MeV 

MeV 

Set 1 

340 

5.595 

2.212 

620.9 

244.3 

Set 2 

400 

5.833 

2.440 

587.9 

240.9 


We should observe that in the eB = 0 limit, ShFR and MFIR regularization schemes 
result in the same mass. On the other hand, the masses for SmFR are about 10 MeV larger 
for both sets of parameters. This happens because the regulator function has a non-zero tail 
for large momentum with respect to the ShFR case, which causes the contribution from the 
vacuum to be larger, thus giving rise to a somewhat larger dressed mass. 


III. NUMERICAL RESULTS 

A. Comparison between regularization schemes 

To motivate the introduction of the MFIR scheme, in this subsection we will compare the 
resulting predictions with the ones of the ShFR and SmFR. For dehniteness, we will consider 
set 1 and take H/G = 0.75, value that follows from various effective models of quark-quark 
interactions ^ . In addition, some general features of the results will be described. 

In Fig. [H we plot the results for M and A as a function of /i, for all three regularizations. 
At a hxed magnetic held, these would appear to present similar behaviors. Two distinct types 
of phases exist in all cases: on the one hand, for low chemical potential, chiral symmetry is 
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broken and snperconducting effects are absent. Since M > /i, all matter terms are zero and 
the dressed mass is independent of the chemical potential. This phase, to be denoted as B 


phase as in previous studies 


[l5, 22- 


25l | , always exists for a low enough chemical potential. If 


chemical potential is increased, on the other hand, there will be a hrst order phase transition 
(whose critical chemical potential pc depends on the regularization) to a phase where A is 
non vanishing and the dressed mass is small. It is considered to be a restored symmetry 
phase, even though exact symmetry restoration occurs only if rric = 0. Since in our case the 
current quark mass is finite but small, restoration is only approximate. For higher chemical 
potentials, further transitions may appear within the A 7 ^ 0 region, causing it to have a 
substructure consisting of several phases. One of such transitions can be seen in the upper 
left panel of Fig. [H signalled by a small kink in the mass at Pc = 334.9 MeV. To understand 
the origin of these we must recall that, for the quark species with |g| = 1/2 and 1 , the 
dispersion relations acquire a Landau level structure due to the magnetic field, as seen in 
Eq.(l7]). Let us consider first the case |g| = 1, where in the matter contribution to the 
thermodynamic potential a sum over these Landau levels up to k^ax has to be performed. 
Here, kmax is determined by the chemical potential, mass and magnetic field. Namely, 
kmax = Floor[(/r^ — M^)/(2ei?)]. Following the notation of previous works, a chirally restored 
phase where the LL’s are populated up to a given k will be referred to as an phase, even 
though we bear in mind that this phase is qualitatively different in that there is a hnite 
diquark gap now. The kink in Fig. [T]then corresponds to a transition in which the highest 
LL that is populated (given by kmax) changes in one unit. In what follows, we will refer to 
these as “van Alphen-de Haas (vAdH) transitions”. It is also clarifying to mention that the 
dressed mass vanishes in the chiral case, so the vAdH transitions are actually signalled by 
discontinuities in the density, and that for the LL they are simply given by the relation 
p = ^j2k\q\B. In the non-chiral case there will be a small departure from this relation 
originating from the finite mass in the restored phase. It is important to stress that out of 
the quark species with different values of \q\, the only one that produces vAdH transitions 
is |g| = 1. In fact, |g| = 0 quarks are decoupled from the magnetic field, so their dispersion 
relation is the same as in the zero magnetic held case. Moreover, quarks with |g| = 1/2 
have an altogether different behavior. The coupling of this quark species to the A removes 
the theta functions from the sum, much in the same way that a theta function becomes a 
Fermi-Dirac distribution when temperature is introduced. Hence, there is no cut off in the 
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sum over LLs, which means that when A is hnite there is non-zero density for all levels. 

Now, even though the order parameters have similar behaviors at a hxed eB as a function 
of chemical potential, the three regularizations exhibit important qualitative differences 
along the magnetic held axis, as can be seen in Figs. [2] and |3l Regularization schemes 
ShFR and SmFR exhibit non-physical oscillations, whose origin lies in the magnetic held 
dependence of the regularization in the vacuum term, that causes the contribution of a 
given LL to be larger for lower magnetic helds. This can be most clearly appreciated in 
Fig. m which displays the behavior of M as a function of eB for p = 0. Here, the only 
contribution to the thermodynamic potential comes from the regularized vacuum. In the 
ShFR, which is the most extreme case, the only LLs participating in the sum are those for 
which > 2k\q\B. Hence, depending on the magnetic held, more or less terms appear and 
each time the relation is satished for a given fc, there will be a discontinuity in the derivative 
of the thermodynamic potential. This singularity, hence, does not correspond to a phase 
transition. The soft regulator, which could be regarded as a way to handle this problem and 
remove sharp oscillations, still contains this pathology, because in this case the contribution 
of a given LL also depends on the magnetic held through the Fermi-type regulator function. 
So, even though the smooth integrals partly conceal this problem, the oscillations are still 
present and in Fig. [2] we can actually see that for ShFR and SmFR they are in phase. On 


the other hand, in the MFIR scheme the mass increases steadily wit 


displaying the usual “magnetic catalysis ehect” as in e.g. Refs. [14, 


22 


1 the magnetic held 


M, 


26|. 


The behavior of the vAdH transitions can also be appreciated in Fig. [3l where we plot 
M and A as functions of eB for /i = 400 MeV. In the range eB = 0.01 — 0.1 GeV^ there is a 
set of peaked discontinuities, each of which corresponds to a vAdH transition. It should be 
emphasized that these are physical transitions, as opposed to the discontinuities previously 
discussed, since they correspond to values of eB at which the quark density for a given LL 
changes from zero to a hnite number. Note that in the ShFR scheme they are harder to 
see because the non-physical oscillations originating from the vacuum contribution are of 
the same order of magnitude. In the MFIR scheme we can also observe that within each 
phase, for a given hnite k, the mass tends to decrease when eB increases, but it shows a 
small increase before the next jump. In the k = 0 phase, the mass decreases steadily with 
magnetic held as well. In the chirally restored phase, the superconducting gap will not 
vanish, and its behavior as a function of magnetic held is nontrivial, as shown in Fig. |3l 
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Also, note that A is approximately constant in the range of eB < 0.12 GeV^, with small 
oscillations resnlting from its conpling to the mass, and then presents a well-shaped curve. 

The phase diagrams in the eB—fi plane as obtained using the three different regularization 
schemes are displayed in Fig. 01 In the case of the ShFR and SmFR, we see that the 
oscillations in the order parameters induce oscillations in the critical chemical potential. 
These are small for low magnetic helds, but become larger in the intermediate eB range and 
once again make the phase diagram hard to interpret. As a result, we can also conclude 
that comparing the three regularizations for a given eB, as was done in Fig. [H is actually 
misleading, since the oscillating behaviors in the order parameters and the critical chemical 
potential in ShFR and SmFR can cause the results to look quite different even for magnetic 
helds which are slightly different. In the MFIR scheme, we note that the critical chemical 
potential is approximately independent of eB for values below 0.07 GeV^, then it decreases 
until it reaches a minimum near eB = 0.2 GeV^ and after this value, it increases back 
again, giving rise to the usual well-shaped curve related to the “inverse magnetic catalysis 


effect’' 


27|. Due to the regulation artifacts, this feature is much less evident in the ShFR 


and SmFR. Goncerning the vAdH transitions, which are the near vertical lines, we note 
that they are almost equal for the three prescriptions. In Fig. 0]they actually correspond 
to the MFIR case, but we make the observation that near the chiral restoration transition 
small deviations exist, which occur because the value of M that produces the deviation with 
respect to \/2keB is different in each scheme. 


B. MFIR results for different model parameters 

In this subsection we further analyze the results obtained within the MFIR scheme, 
paying particular attention to their dependence on the model parameters. Results for the 
two parametrizations introduced in Sec.III.G will be given. Moreover, in the previous section 
only H/G = 0.75 was considered. However, given that the value of this ratio is subject to 
certain degree of uncertainty, it is worthwhile to explore the consequences of varying it 
within a reasonable range. Thus, in what follows, the representative values H/G = 0.5, 0.75 
and 1 will be considered. A few comments on how the model results change for H/G < 0.5 
will be also made. Note that values H/G > 1 are quite unlikely to be realized in QGD. 

Let us start by analyzing the behavior of the order parameters as a function of /i for 
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given values of magnetic field. We will concentrate on the results obtained with set 1, since 
it exhibits a more complex phase structure. Set 2 will be addressed further on. As was seen 
in the previous section for H/G = 0.75, the system is in the B phase for low /i, where the 
dressed mass is large. On the other hand, it is in one of the possible A-type phases for a high 
enough fi value where the dressed mass is small. However, if the coupling ratio is changed, 
other phases may appear for low eB values and intermediate chemical potentials. This can 
be seen in the left panels of Fig. [5l where eB = 0.04 GeV^. For H/G = 0.5, there is at 
/ic = 338.2 MeV a weak hrst order transition from vacuum to a phase where the mass is 
slightly lower and also a slowly decreasing function of /r. Quark density is hnite for the |^| = 0 


and |q| = 


Refs. 


15 


22 


quarks (only the lowest LL being occupied for the latter species). Following 


25| this phase will be denoted as a C-type phase although, as in the case of the 


A-type phases, here the superconducting gap is non zero. Actually, it happens to be very 
small, remaining always under 1 MeV. Hence, it is not visible in this scale. If /x is further 
increased, we hnd another hrst order transition to an A-type phase, at 346.1 MeV. However, 
we note that A is in the range 25 — 30 MeV, which is a relatively small value compared to the 
resulting ones from higher coupling ratios in A-type phases. As the coupling constant ratio 
is increased from the value H/G = 0.5, the upper phase transition displaces downwards. 
This causes the C-type phase to shrink until it eventually disappears around FT/G ~ 0.65, 
so that a single phase transition remains connecting the B phase to the A-type phases. As 
the coupling ratio is further increased, this phase transition continues to move downwards 
and for H/G = 0.75, the transition occurs at /ic = 333.6 MeV. For even larger values of 
H/G (~ 0.94), the phase transition splits into two once again, so that for H/G = 1 there is 
an intermediate phase, which will be referred to as a D phase. This phase is qualitatively 
different from the one found in H/G = 0.5. To begin with, the transition from vacuum to 
this phase is second order. The dressed mass is still large but the superconducting gap is 
hnite and actually increases sharply with /i. Since both condensates are appreciably large, 
this is usually referred to as a “mixed phase” 28|], even though other meanings exist in the 
literature for this term [29l]. There is no quark population for the |g| = 0 and |g| = 1 species, 
but the hniteness of A induces a non-zero density for \q\ = 1/2 quarks. The transition 
leading to the A-type phases is hrst order as in the previous cases and A ~ 175 MeV. 
Within it, there is another transition for /i = 308.5 MeV, which is actually a vAdH being 
traversed vertically. The behavior is much simpler for eB = 0.3 GeV^, as can be seen in 
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the right panels of Fig. [5l For all values of the coupling ratio, we only see the B and Aq 
phases. Even though transitions to higher phases will appear for much higher chemical 
potentials, it is generally seen that the phase structure is simpler for eB > 0.15 GeV^, so that 


a single transition connecting the B phase to the Aq phase is seen in the range of /i which is 


of interest (see 


23| for a detailed discussion on how this occurs in NJL with magnetic field 


and without diquark pairing). For both magnetic helds, we see that increasing H/G always 
produces larger values for A and reduces the fic necessary to achieve the superconducting 
phase. We should also note that for smaller magnetic helds, A grows faster as a function of /i 
than for eB = 0.3 GeV^, where it seems to be almost constant. As a matter of fact, we also 
checked that for magnetic helds larger than a value around 0.4 GeV^, the superconducting 
gap actually decreases (yet only slightly) with chemical potential. 


A deeper understanding of the behavior of the phases can be obtained from the phase 
diagrams in the eB — /i plane (Fig. |6]), where our three coupling ratios and both parameter 
sets are considered. The hrst transition encountered if the phase diagram is traversed in the 
direction of increasing fi will be referred to as the “main transition”. It has approximately 
the same shape in all displayed phase diagrams and it connects the B phase to populated 
phases in general. In set 2, the fact that Mq adjusts to a higher value causes the main 
transition to occur at a higher chemical potential. On the other hand, as H/G increases, 
the main transition is displaced downwards in its entirety, and the depth of the “inverse 
magnetic catalysis well” diminishes. In set 1, we can see for H/G = 0.5 (top left panel) 
that both transitions enclosing the G-type phase are constant in a large magnetic held range 
so that it extends in an approximately horizontal band up to eB ~ 0.11 GeV^, where it is 
bounded by a crossover type transition, signalled by the peak of the chiral susceptibility. The 
crossover leads to an Aq phase and, as expected, the mass drops sharply and A increases. 


The different possible criteria to dehne this kind of transitions were discussed in detail in 


23|. 


On the other hand, if magnetic held is decreased down to zero, the two horizontal transitions 
continue to exist (however, the lower transition becomes second order when eB = 0), which 
means that a phase with large mass and quark population exists in the NJL model without 
magnetic held. The detailed behavior of the G-type region is rather involved in the limiting 
case in which eB tends to zero, and will not be discussed any further. It is important to point 
out that the existence of G-type phases is not a consequence of the diquark pairing channels 
either. As a matter of fact, up to H/G = 0.5 superconducting ehects are still small, and the 
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phase diagram remains basically unchanged. In fact, the two roughly horizontal transition 
lines are seen to remain almost unmodihed if the coupling ratio is swept between these two 
parameters. We observe that even though this C-type region consists of a single phase where 


only the lowest LL is populated (hence it is a Cq phase), the results from [23|, |2^ suggest 


that it may become more complex if vector interactions are introduced, or for parameter sets 
leading to a lower Mq value. In that case, it could be possible that phases with higher fc’s 
appear between the two horizontal transitions. In the A-type phases, higher LL population 
is allowed because the mass is lower than in the C-type phases. vAdH transitions separate 
phases with different number of LL populated, and in the presented diagrams, LL’s up to 
k = 7 are occupied by |g| = 1 quarks. On the other hand, the phase diagram is simpler in 
top right panel), where it is seen that there is no intermediate phase for H/G = 0.5. 


set 2 
In 


231] we showed that changing the parameter set so that Mq increases always produces a 


simpler phase structure where there are only a main transition and vAdH transitions. As 
has been said, the upper transition found in set 1 moves downwards when H/G is increased, 
so the Co phase shrinks and hnally disappears, leaving a rather simple phase structure for 
H/G = 0.75 (middle-left) panel. Note that for this value of the coupling ratio the phase 
diagrams for both sets are qualitatively similar. Finally, the phase diagrams for H/G = 1 
are shown in the bottom panels, where the D phase is seen to exist for both parameter sets. 
Once again, we observe that the existence of this phase is a consequence of the diquark 
pairing alone and hence it already exists for eB = 0. It is interesting to note, however, 
that it continues to exist for hnite magnetic helds. The two transitions delimiting it are 
roughly horizontal for low eB. For higher values, the transitions lines move closer together 
and hnally intersect near 0.1 GeV^ causing the D phase to disappear. At low magnetic held, 
the D phase exists for a narrow /i range of at most 5 MeV for either set, but larger values 
oi H/G would cause the phase to extend farther in the chemical potential direction. 


In Fig. [71 we present the behavior of the order parameters as a function of magnetic 
held, for all considered values oi H/G and representative values of fi. Once again, we 
display the results for set 1 only, since no qualitatively diherent behaviors arise for set 2. 
It was already seen in Sec.III.A that in the B phase, mass increases with magnetic held 
and that A = 0. This result, which is a manifestation of the magnetic catalysis ehect, is 
exclusively seen in vacuum. Since A = 0 in this phase, the behavior is independent of H/G. 
When quark population is hnite, several possibilities arise depending on the magnetic held. 
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parameter set and coupling constant. For high chemical potentials (/r = 360 MeV), where 
chiral symmetry is restored for all values of H/G, the mass exhibits a series of peak-like 
discontinuities. They correspond to the already discussed vAdH transitions, where the k^ax 
corresponding to the |g| = 1 quark changes in one unit. Due to the coupling to the mass, A 
also oscillates. These oscillations become particularly large for intermediate magnetic held, 
around eB = 0.1 GeV^. We note that on increasing H/G, M is shifted downwards and A 
is shifted upwards. Also, all oscillations are smoothed and become negligible compared to 
the scale of the order parameters. 

For H/G = 0.5, fi = 340 MeV corresponds to the intermediate C-type phase where quark 
population is hnite and A is small. We note that for eB near zero, the mass is close to the 
vacuum value, but gradually decreases as eB is increased. This feature is a manifestation of 
the anticatalysis effect. On the other hand, A is very small and increases only slightly, always 
remaining under 1 MeV. The jumps in both M and A near eB = 0.1 GeV^, correspond to 
a hrst order transition to an Ai phase. For larger magnetic helds, the mass decreases 
steadily until it hnally returns to the B phase, which corresponds to the high magnetic held 
branch of the main transition. In the lower two panels of Fig. [TJ we also see the behavior 
of M and A in the mixed phase for H/G = 1. Here, we can see that the mass is almost 
independent of magnetic held in the D phase (slightly decreasing), while the superconducting 
gap increases from 25 MeV to 50 MeV. A few comments should be made regarding the vAdH 
transitions. These transitions are a characteristic of the chirally restored region of any NJL 
model with magnetic held. In particular, in the model without diquark pairing, there will be 
vAdH transitions for up quarks and down quarks. In the rotated base, quarks with charges 
|g| = 1 and \q\ = 1/2 couple to the magnetic held, of which only the former yield hrst order 
transitions. However, there is a vAdH-like cross over transition for H/G = 0.5 and 0.75 that 
resembles a \q\ = 1/2 vAdH in the following sense. If we were to set H = 0, then A = 0, 
and the integrals in the |g| = 1/2 matter terms would transform to Heaviside functions, 
cutting oh the LL sum at a maximum k given by the relation iH — (M^ -|- keB) and 
therefore producing vAdH transitions. If if/G is made hnite but small, these discontinuities 
quickly disappear, but a remnant of these transitions is still present because the associated 
quark number susceptibility still exhibits peaks. However, they are smeared out as H/G is 
increased and for H/G = 0.5 there is no trace left of these transitions, except for the one 
corresponding to the LL passage from 0 to 1. In set 1, we can see that this transition is 
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first order and has an end point after which it becomes a cross over, as can be seen in the 
corresponding panel in Fig. [6l This crossover is still present for H/G = 0.75 and hnally 
disappears for H/G = 1. We make the observation that since these are rather weak, we do 
not consider them to separate distinct phases. 


IV. SUMMARY AND CONCLUSIONS 


In the present work we explored the effects of magnetic held on cold color superconduct¬ 
ing quark matter in the framework of the NJL-type model, using a regularization scheme 
in which the contributions which are explicitly dependent on the magnetic held turn out to 
be hnite and, thus, do not require to be regularized. Such a “magnetic held independent 
regularization” (MFIR) scheme can be considered an extension of the method described 
in e.g. Ref.Q to the case in which color pairing interactions are present. We compared 
the correspondin g re sults with those obtained through the regularization methods used in 


previous works 


13|. In those works, a regulator function was introduced in order to sepa¬ 


rately render the contribution of each Landau level hnite. This, however, might lead to the 
appearance of non-physical oscillations even when rather soft regulator functions are used, 
while within the MFIR scheme these oscillations are completely removed and, thus, results 
are easier to interpret. Note that in the cases we investigated the smooth regulator case 
exhibits barely no diherence with respect to the MFIR scheme for eB up to 0.1 GeV^, and 
the above mentioned oscillations start to become relevant from this value on. In respect of 
this, it should be borne in mind that some estimates}^ indicate that the magnetic helds at 
the center of magnetars can be as large as eB ~ 0.6 GeV^. The remaining of this paper was 
devoted to investigate, for the MFIR regularization scheme, the model parameter depen¬ 
dence of the behavior of magnetized cold superconducting quark matter. We considered two 
parameter sets that adjust to acceptable values of the dressed masses, and that were already 
cnown to generate qualitatively different phase diagrams in the non-superconducting case 


23l |. Moreover, three representative values of the coupling constant ratio H/G were consid¬ 


ered. We found that up to H/G = 0.5 superconducting effects are still small, and the phase 
diagram remains basically unchanged. In particular, for set 1 there is an intermediate phase, 
in which quark population exists but chiral symmetry is still strongly broken. The diquark 
gap in this phase is hnite but extremely small. It is connected to other phases by hrst order 
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transitions and a chiral crossover. As the conpling ratio is increased, the two transitions 
snrronnding this phase move close to each other, cansing it to disappear at H/G ~ 0.65. A 
relatively simple phase diagram hence exists for a narrow range aronnd the standard valne 


H/G =0.75. If the conpling ratio is increased beyond 0.94, a mixed phase is present in 


both parameter sets, where both condensates are large. This phase already exists for zero 
magnetic held and we showed that according to this model it extends in the magnetic held 
direction for a relatively large range of ei?, disappearing for eB O.lGeV^. It is bonnded 
from below by a second order transition. The natnre of the vAdH transitions was also 
discussed. We explained that in the rotated base these exist only for the |g| = 1 quark, 
and that even though the |g| = 1/2 quark couples to the magnetic held, its coupling to the 
superconducting gap smears the vAdH transitions out, leaving only a crossover remnant for 
the transition between the two lowest k values. Finally, we found that the inverse catalysis 
phenomenon is observed in all the cases, although less pronounced as H/G increases. 

Throughout this work we have concentrated on the analysis of the impact of the novel 
regularization scheme on the model predictions for the ehects of magnetic held on symmetric 
two havor cold color superconducting quark matter. It is clear that, for applications in 
the physics of compact stars, the neutrality and /5-equilibrium conditions should be taken 
into account. Moreover, to address the ehect on the description of the CFL-type phases 
strangeness degrees of freedom have to be included. We expect to report on these topics in 
forthcoming publications. 

APPENDIX 

In this Appendix we present some details concerning the derivation of Eq. ljTTll . We start 
by considering Eq.(I7j). As a hrst step towards the regularization of this expression we sum 
and subtract the contribution for vanishing chemical potential. In this way we get 


P\q\ = l/2 — Si + S 2 


(13) 


where 



(14) 


and 



(15) 
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where 


(+ M2 + s /i)2 + A2 _ Vz + M2 + A2 


s=±l 


( 16 ) 


It is clear that both Si and S'2 are divergent. However Si has the standard form analyzed in 
Ref.[^, with m2 —)■ A2 + M"^. Thns, following the steps described in that reference we get 

{eBf 


Si = ^ [ dp V^p2 A2 + 


TT^ 


27r2 


c>( 1 \ I 

^ (-1,2/) + J 


(17) 


where y = (M2 + ^^)/eB and, as in Ref. Q, a 3D sharp cntoff has been introdnced to 
regnlarize the eR-independent contribntion. 

To regnlarize S'2 we add and snbtract the corresponding contribntion in the absence of 
magnetic held. Namely, 


^2=4 


d^p 

{2Tlf 


F{p^) + R 


where 


eB 


k=Q 


djh 

2 tt 




d?p 

Wy 


F{p^ 


(18) 


(19) 


We prove in what follows that R is hnite. For this pnrpose we start by rewriting its second 
term introdncing cylindrical coordinates. Then, once the angnlar integral is performed we 
get 


d?p 


dp. 


dp 


eB 


dpz 


FiP ) = ^ if p nvi + PY = -r tP / dx F{pi + eBx) (20) 


27r 


27r 


dvr 


27r 


where the change of variables p = \/eBx has been nsed in the last step. Replacing this 
expression in Eq. flT^ we get 


eB 


dp. 


eB 


dpz 




27r 


2 ti 


TT 


27r 


X; F(pI + klB) 


dx F{pI + eBx) 


k=l 


( 21 ) 


For convenience we introdnce at this stage f{z) = F{eBz)/\/eB and perform the change of 
variables p = Pzj^/eB. Then we get 


R = 


{eBf 

27r2 


dp f{pY + 


(e~R)^ 


TT^ 


dp 


5^/(/ + ^) 


dx f{p^ + x) 


k=l 


( 22 ) 


where the fact that the integrands are even fnnctions of p has been nsed. To proceed we 
notice that f{z) is bonnded in the interval 0 < z < 00. Thns, to prove that the momentnm 
integrals are convergent it is enongh to verify that the corresponding integrands vanish 
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sufficiently fast as p ^ oo. For this purpose it is convenient to consider the expansion of 
f{z) for ^ >> /i^/ei?; (M^ + A‘^)/eB. We get 

M = Z 7^72 (23) 

m=l 


where the values of Cm for m = 1, 2 are 

fA _ 3M2 +A2' 
eB 2 eB 


Ci = 


Wr^ 


c,= 


A^/i" 


(24) 


It is clear from Eq.(|23]) that for large p the integrand in the first term of R goes as 1/p^ 
and, thus, the corresponding integral is convergent. To prove the convergence of the second 
term one can proceed as follows. Let us consider N » pP'/eB] (M^ + A?)/eB such that 
ioi z > N we can replace f(z) by the expansion Eq. (l23ll . As an example, we note that 
for the range of values of p and eB considered in this work, /eB and (M^ + A^)/ei? are 
always smaller than 25. Accordingly, we have checked that if A >> 25 (e.g. N=250) the 
use of the hrst few terms in the expansion Eq. fl23l) leads to values which are in excellent 
agreement with those obtained from the full expression. Then, denoting by I the integrand 
in the second term of Eq. fl2^ , we have 

' N „Ar 


I 


V f{p^ + k) - dx f{p^ + x) 
L.-1 Jo 


+ Cm 

m=l 


oo 


E 

_k=N+l 


1 

(p2 + fc)m+l/2 


dx 


In {p^ + a;)’”+F2 


(25) 


Next we note that, since m > 1, the sum over k in the second term can be written in terms of 
the Hurwitz zeta function ^ and that the integral over x is convergent and can be explicitly 
performed. Consequently we get 

' N ^N 


I 


V f{p^ + k)- dx f{p^ + x) 
1.-1 Jo 


m=l 


^ (m + 1/2,1 + + a) - 


777 ,- 1/2 (p2 -|- 


(26) 


Considering values of >> A and making again use of Eq. (|23|l it is easy to see that the 
two terms in the hrst bracket of Eq. fl26ll go as A/p^ which implies that their contributions 
to R are hnite. What remains is to show that the leading contribution to the second bracket 
also goes faster than 1/p. This is not so obvious since for m = 1 the second term does go 
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as 1 /p. Actually, as shown in what follows we need this term to go like that in order to 
cancel a similar term arising from the hrst term. In fact this is the main reason why the 
substruction scheme proposed in this work was introduced. To show how this cancellation 
occurs we consider the expansion of the Hurwitz zeta function. For a; >> 1 and m > 1 one 
has 

-(m-l/2) -(m+1/2) 

«m + 1/2,1 + x) =-—---+ O (27) 

m — 1/2 2 ' 

Using this expansion it is easy to see that for p^ » N we have 


i{m + 1/2, l+p'^ + N) - 


+ O (28) 


m-l/2 (p2 + iV)”^-V2 ' 2 p 2 ™+i 

This implies that the leading contribution comes from the m = 1 term and it goes as 1/p^ 
ensuring the convergence of the corresponding integral. 

Having proved that R is indeed hnite only the hrst term in S 2 requires regularization. 
Using a 3D sharp cutoff as before we get 
2 




TT^ 


dp p 


Ei + EX- ‘l-^z + M2 + A2 




27r2 


^00 noo 

/ dp ^ak /(p^ + k)-2 dx f{p^ + x) 
Jo Jo 


(29) 


where, as dehned in the main text, E^ = y {\/p^ + ± ^)2 -|- Finally, replacing 

Eqs. (11711291) in Eq. flT^ we immediately get Eq. dTT]) . 
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FIG. 1. Left panels: M vs /x (upper) and A vs /x (bottom), for eB = 0.052 GeV^. Right panels: 
M vs fi (upper) and A vs /x (bottom), for eB = 0.30 GeV^. Set 1 was used, H/G = 0.75 and all 
regularization schemes considered. 
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FIG. 2. M vs eB in the B phase, for the three regularization schemes considered, set 1 and 
H/G = 0.75. Note that A = 0 in this phase. 
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H/G=0.75 - n = 400 MeV 
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FIG. 3. M and A vs eB in the A-type region, for the three regularization schemes considered, 
set 1 and H/G = 0.75. 
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FIG. 4. Phase diagram in the eB vs n plane, for the three regularization schemes considered, set 1 
and H/G = 0.75. All transitions seen in this diagram are first order. 
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FIG. 5. M and A vs /x within the MFIR scheme for set 1, different H/G ratios and two values of 
eB. Full, dashed and dotted lines correspond to H/G = 0.5,0.75 and 1 respectively. 
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FIG. 6. Phase diagrams within the MFIR scheme for different H/G ratios, for both set 1 (left 
panels) and set 2 (right panels). Full black lines correspond to first order transitions and dotted 
black lines to second order transitions. The dash-dotted red lines represent chiral susceptibility 
crossovers. Blue dashed lines represent quark number susceptibility crossovers which are too weak 
to be actually considered to separate distinct phases. Note that the sets are plotted in different 
intervals of the /r axis, even though the scale is the same. 
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FIG. 7. M and A vs eB within the MFIR scheme for the three values oi H/G considered and 
several representative values of /r, set 1. 

































































